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The discovery of the high temperature superconductivity in copper oxide materials has 
initiated an enormous theoretical effort on quantum antiferromagnets in two- dimensions 
(2D) 0-0. Although experiments have not yet isolated the key elements of the electronic 
structure necessary for a global understanding of the physical properties of copper oxide 
materials, a significant body of reliable and reproducible data has been amassed by using 
many probes HHH. The normal-state properties of these materials exhibit a number of 
anomalous properties in sense that they do not fit in the conventional Fermi-liquid theory 
g. When undoped, with one hole per copper site, these materials are antiferromagnetic 
Mott insulators ||. When doped with sufficient carriers, they are superconductors with 
no the antiferromagnetic long-range order ||. There is strong evidence from transport and 
neutron scattering that no sharp phase transition occurs, and the strong correlations are very 
important to the electronic structure ||. The angle-resolved photoemission spectroscopy is 
the definitive way to determine the electron dispersion relation for any materials. In copper 
oxide materials, the angle-resolved photoemission experiments || have produced interesting 
data that introduce important constraints on theories. It has been shown [|2]|| that the hole- 
doped copper oxide materials exhibit universal properties likely induced by the behavior of 
carriers in their common copper oxide sheets. It is believed [TJ that the physics of these 
materials may be effective described by a 2D, large U, single-band Hubbard model. In the 
large U limit, the Hubbard model is transfered into t-J model acting on the space with 
no doubly occupied sites. Furthermore, Zhang and Rice |J derived the t-J model from a 
multiband Hubbard model described 2D copper oxide planes. 

It is clearly very important to establish appropriate formalism for the problem and to 
show that this leads to behavior similar to that seen in experiments. In order to account for 
real experiments under the t-J model, the crucial requirement is to impose the electron local 
constraint 0. The local nature of the constraint is of prime important, and its violation may 
lead to some unphysical results || . Recently a fermion-spin theory based on the charge-spin 
separation is proposed || to incorporate this constraint. In this approach, the electron on- 
site local constraint for single occupancy is satisfied even in the mean-field approximation 
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(MFA). In the framework of the fermion-spin theory, the ground-state properties, such as, 
the ground-state energy, ground-state kinetic energy, phase separation, specific heat data, 
and doping dependence of the antiferromagnetic long-range order, are discussed |p|-|TTH and 



the results are in qualitative agreement with experiments and numerical simulations. In 
this paper, we develope a mean-field theory in optimal doping regime within the fermion- 
spin theory to study the photoemission spectrum, the electron dispersion, and the electron 
density of states, which is useful for understanding of the electronic structure of the copper 
oxide materials. 

We start from the t-J model which describes the electrons moving on a planar square 
lattice, 

H = -tJ2 clc ja + h.c. -fij2 c l c ^ + J E 8* ■ s jt (i) 

(ij)a icr (ij) 

where C\ a (Ci a ) are the electron creation (annihilation) operators, Sj = CjaCi/2 are spin 
operators with a = (a x , a y , a z ) as Pauli matrices, and fi is the chemical potential. The 
summation (ij) is carried over nearest nonrepeated bonds. The t-J Hamiltonian (1) is sup- 
plemented by the on-site local constraint, S<j Cj^C^ — 1> ^- e -> there be no doubly occupied 
sites. With the help of the fermion-spin transformation || 

Ci\ = h\S~, Cii = h\S+, (2) 

where the spinless fermion operator hi keeping track of the charge (holon) while the pseu- 
dospin operator Si keeping track of the spin (spinon), the t-J model (1) can be rewritten as 

1 

H = -tJ2hih}(SfSr + SrS+) + h.c. 

(ij) 

-fiJ2h\hi + jJ2(hM)(Si ■ SjXtyiJ), (3) 

* (ij) 

where and S~ are pseudospin raising and lowering operators, respectively. It is shown 
0] that the constrained electron operator can be mapped exactly onto the fermion-spin 
transformation defined with an additional projection operator. However, this projection 



operator is cumbersome to handle for the actual calculation possible in 2D, we have dropped 
it in Eq. (3). It has been shown in Ref. || that such treatment leads the errors of the order 
5 in counting the number of spin states, which is negligible for small doping 5. Within the 
MFA, the t-J model (3) can be decoupled as, 

H = H t + Hj- 8Nt X <p, (4a) 
H t = 2 X t £ h\ +v h -nY, h\h h (4b) 

Hj = -^Jeff^Y^t^i+r, + S i S t+ v ) + Jeff Y ( 4c ) 
i,r] i,'<] 

with r] = ±x,±y, N is the number of sites, and J e jf = J[(l — 5) 2 — <p 2 }. The nearest- 
neighbor spin bond-order amplitude x and holon particle-hole parameter (ft are defined as x — 
(S^Sf +v ) and = (h\h i+v ), respectively, where the site subscripts of the order parameters 
X and 4> have been dropped since the system is translation invariant. In this mean-field 
level, the spinon part is described by an anisotropic Heisenberg model with the anisotropic 
parameter is given by, 

£ = J eff + 2t0 ^ 
Jeff 

The quantum spin operators obey the Pauli spin algebra, and this problem can be dis- 



cussed in terms of the two-time spin Green's function within the Tyablikov scheme [12]. In 



this case, the one-particle spinon and holon two-time Green's functions are defined as, 

D(i-j,t-t) = -id(t-t')([St(t),Sr(t')}) = ((St(t);S~(t'))), (6a) 

D,(i -j,t- f) = -i6(t - t')([S*(t), S*(t')\) = ((S*(t); S*(t'))), (6b) 

and 

g(i- j,r-r') = -i0(t-O<[M*).*}(O]> = ({hi(t);h](t'))), (7) 

respectively, where (...) is an average over the ensemble. Because the spinon system is an 
anisotropic, we have defined the two spinon Green's function D(i—j, t—t') and D z (i—j, t—t') 
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to describe the spinon propagations. The time-Fourier transform of the two-time Green's 
function satisfies the equation, 

W ((A;5)) = i([AB]> + (([Aff];B)>, (8) 

therefore the Green's functions can be obtained by applying the Tyablikov decoupling tech- 
nique, and the correlation functions can be obtained by the spectral representations as 

(B(t')A(t)} = i r ^(^^>W-(^ g >W e -M t -Q. (9) 

J — oo 27T C' — 1 

Recently, we |9HIT| have employed the fermion-spin theory to study the ground-state 



properties of the 2D t-J model and obtained some interesting results. Within the random- 
phase approximation, we have shown that the antiferromagnetic long-range order is 
destroyed by hole doping of the order ~ 5% for the reasonable value of the parameters 
t/J = 5. Then in the following discussions, we only study the systems in the optimal 
doping regime (20% > 5 > 5%), where there is no the antiferromagnetic long-range order, 
i.e., (Sf) = 0. In this case, the basic equations for the mean-field spinon two-time Green's 
function in one-dimension have been discussed in detail by Kondo and Yamaji |[L3|| . Following 



their discussions, we can obtain the mean-field spinon Green's functions of Hj in 2D, 

= A[(*x, + xh»-(<X + 2x,)] (_! 1 ) , (10 .) 

Dz(K W) = ~ 7TT 77T 7TT , (10b) 

2u) z {k) \uj — uj z {k) u + u z (k)J 
where 7k = 4 J2 V e* k '^, and 

u\k) = A 2 (ae( Xz l k - \x){^ - 1) + [aC z + ^(1 - a)](l - e 7 *)) 

+A 2 f-aexjkilk - e) + \e[aC + — (1 - a)](e - 7fc )) , (Ha) 



^(/c) = A + — (1 - a)] - ae X 7fc - ^ae*) (1 - Tfc), (Hb) 

with A = 2ZJ e ff, Z is the number of the nearest neighbor sites, and the order parameters 

Xz — {SfS^), C = ^ X^'^i+^i+r/)) an< l C z = Ysqtf (^i+ri^i+rj') • l n or der not to 



violate the sum rule of the correlation function (S^S^ ) = \ in the case {Sf) = 0, the 
important decoupling parameter a has been introduced as these discussed by Kondo and 



Yamaji ||13|| , which can be regarded as the vertex corrections. At half-filling, the t-J model 
is reduced as the isotropic antiferromagnetic Heisenberg model, where e = 1, Xz — f X? an d 
C z = \C in the rotational symmetrical case, and we obtain D z (h, uj) = \D(k, u), which is 
just these discussed by Shimahara and Takada [14|. 

The mean-field Green's function of H t for holons is very simple, and can be written as 

<?(k,u,) = i -, (12) 

u-(e k - /I) 

where Sk = 2Zxt'jk- With help of the spinon and holon Green's functions and the spectral 
representations of the correlation functions, the order parameters Xi C, Xz, C Z1 0, and 
chemical potential \i can be obtained by the self-consistent equations, 

X = ^ X> A[<2% + X) Z\~ fa + ^' coth^), (13a) 

C = ±X* A[(2eX ' + %) + ^W ^f ). d3b) 

1 = 1 A[(2 % + X)7t - (e X + 2X.)] Mj) 

2 N Y 2u(/c) ( 2 ' v ' 

X.= jyl> 2 ^ w cotht-^-), (13d) 

*=^?4-*^). (1*) 

As we have shown in the previous works [O] that the present MFA self-consistent calculation 
is just the usual self-consistent Hartree-Fock approximation. 

We |H| have performed a numerical calculation for these mean-field self-consistent equa- 
tions. The results for the order parameters at optimal doping regime are very close to our 
previous works HfjlJ based on the 2D Jordan- Wigner approach, the detailed discussions 
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will be given elsewhere [jnj. In this paper, we hope to discuss the electronic structure 
of copper oxide materials, and therefore it need to calculate the electron Green's function 
G(i - j,t - f) = ((Ci a (t);Cj a (t'))). According the fermion -spin transformation (2), the 
electron Green's function is a convolution of the spinon Green's function D(i —j,t — if) and 
holon Green's function g(i — j,t — t'), and can be obtained at the mean-field level as, 

t F ^v) + F ^ ) , (14 ) 

\U - Ul(p) + Ep- k UJ + L0(p) + Ep^k J 

where F 1 (k, p) = n B (u p ) + n F (e p _ fc - (if) , F 2 (k, p) = 1 + n B {u p ) - n F {e p - k - (i) , with n B (uj p ) 
and np{E p -k — A 4 ) are the boson and fermion distribution functions for spinons and ho Ions, 
respectively. From the electron Green's function (14), we obtain the electron spectrum 
function, 

a fir \ ot nn \ o 1 A [( 2 % + X)% - (^ + ^Xz)] w 
A(k, u) = -2ImG(k, u) = 2vr— ^ x 

(F 1 (k,p)5(w-^(p)+£ p _ fe ) + F 2 (k,p)(5(cj + cj(p) + e p _ fc )) . (15) 

In the t-J model, the doubly occupied Hilbert space has been pushed to infinity as Hubbard 
U — > oo and therefore the spectrum function only describes the lower Hubbard band. Our 
mean-field result of the spectral functions at the doping 5 = 0.12 for the parameter t/J = 2.5 
is shown in Fig. 1 (solid line). For comparison, the exact diagonalization and quantum Monte 



Carlo result at 5 = 0.12 for t/J — 2.5 obtained by Moreo et al. [161 is also shown in Fig. 1 
(dashed line). Although the particular details of the spectral function and dispersion may 
differ from compound to compound, some qualitative features seem common to all copper 
oxide materials. Hence a quantitative comparison between theory and experiment is still 
early, but the qualitative tendency of the spectral function and dispersion in an adequate 
theoretical description should be consistent with experiments and numerical simulations. 
In the present mean-field theory, the most important feature is that the intensity peaks 



is qualitative consistent with the numerical simulation [ lEfl . The low energy peak is well 
defined at all momenta, and the positions of the dominant peaks in A(k, u) as a function of 



momentum are shown in Fig. 2, which is also in qualitative agreement with the numerical 



simulation [IB] and the experimental result 



Now we consider the electron density of states, which is defined as 

p(w) = iE^ w )- ( 16 ) 

k 

The numerical analysis of the electron density of states in the t-J model as a function of 
doping has been done by many authors ||17| . On the other hand, oxygen x-ray absorption 



spectra measured |L8j on La2- x Sr x CuO^ may be interpreted in terms of a picture in which 
hole doping introduces carriers into the lower hand. Our mean-field result at doping 5 = 0.12 
and 5 = 0.06 for the parameter t/J = 2.5 is shown in Fig. 3. We find that the chemical 
potential fi moves from nearly zero at small doping 5 = 0.06 to the top edge of the lower 



Hubbard band, which also is qualitative consistent with the numerical simulation [XT] and 
experimental result [|l^ . 



In summary, we have developed a mean-field theory which satisfying the electron on-site 
local constraint in the relevant regime of density for the high temperature superconductors, 
namely in the vicinity of optimal doping within the fermion-spin theory. Within this mean- 
field theory, we have study the electron spectral function, the electron dispersion, and the 
electron density of states of copper oxide materials, the results are qualitative consistent 
with the numerical simulations and experiments. 

At zero doping, the t-J model reduces to the antiferromagnetic Heisenberg model, which 



has a local SU(2) symmetry in the fermion representation [ 19[ . This symmetry does imply 
that the spinon particle (hole) state with spin up is the same state as a spinon hole (particle) 
with spin down. In the conventional slave-boson theory |p0[| , the SU(2) is broken to U(l) 
upon doping, and this U(l) gauge degree of freedom is introduced to incorporate with the 
single occupancy local constraint. However, in the fermion-spin theory, the SU(2) symmetry 
is also broken upon doping. Moreover, since the local constraint is satisfied exactly even in 
the MFA, the extra gauge degree of freedom occurring in the slave-particle approach does 
not appear in the here, which is consistent with these discussions in Ref. ||21|| . 
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Final we also note that recently Wen and Lee [22] have developed a slave-boson theory 
for the t-J model at underdoped regime which preserves SU(2) symmetry, they argued that 
spin gap phase at the underdoped regime can be understood as the staggered flux phase. 
However, the area of the Fermi surface produced by this SU(2) mean-field theory is larger 
than the predicted by the Luttinger theorem which reveals a drawback of the SU(2) mean- 
field theory at the optimal doping, and therefore the U(l) mean-field theory is better at the 
optimal doping. 
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FIGURES 

FIG. 1. Spectral function A(k, to) of the 2D t-J model within the fermion-spin mean-field theory 
(solid line) and the exact diagonalization and quantum Monte Carlo methods (dashed line) for the 
parameter t/J = 2.5 

FIG. 2. Position of the the dominant peaks in A(k, to) as a function of momentum. 

FIG. 3. Electron spectral density of the 2D t-J model at the parameter t/J = 2.5 for (a) doping 
5 = 0.12 and (b) S = 0.06. 
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